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A comprehensive model, called ECOMB, is proposed to describe multiparticle produc¬ 
tion by soft interaction. It incorporates the eikonal formalism, parton model, color mutation, 
branching and recombination. The physics is conceptually opposite to the dynamics that 
underlies the fragmentation of a string. The partons are present initially in a hadronic colli¬ 
sion; they form a single, large, color-neutral cluster until color mutation of the quarks leads 
to a fission of the cluster into two color-neutral subclusters. The mutation and branching 
processes continue until only qq pairs are left in each small cluster. The model contains self¬ 
similar dynamics and exhibits scaling behavior in the factorial moments. It can satisfactorily 
reproduce the intermittency data that no other model has been able to fit. 


I. INTRODUCTION 

The study of multiparticle production in low-pr processes has been pursued for over twenty years. Since 
they involve soft interactions, they cannot be treated in perturbative QCD. In the absence of any reliable 
theoretical approach to the problem, many models have been proposed, most of which are represented in 
the review volumes published ten years ago Nearly all of those models have since been shown to be 

inadequate in light of the data on fluctuations and intermittency Indeed, there are very few models 

that have the appropriate dynamical content capable of reproducing the scaling behaviors observed in the 
experiments. Even though perturbative QCD cannot be used, it seems that some kind of branching process 
is needed to generate the property of self-similarity, as the resolution scale is varied. In this paper we propose 
a model that incorporates some aspects of nonperturbative QCD and is capable of generating the features of 
the intermittency data, which are shown in the last hgure of this paper. To our knowledge those data have 
not been fitted by any model that contains some features of the color dynamics. To reproduce those data 
has become the main motivator for this work. 

Our approach embraces many time-honored properties of hadronic collisions. Since hadrons are extended 
objects, the eikonal formalism is at the foundation of our model. Thus it is not difficult for our model to 
possess the virtues of geometrical scaling 1^ and approximate KNO scaling In order to build into the 
model features of chromodynamics, it is necessary to introduce quarks and gluons into the eikonal formalism, 
so the parton model is an essential gateway into the microscopic domain of color interactions. Once we enter 
that domain, we embark on an unconventional journey of studying color mutation of the constituents as a 
dynamical process by which the colors of the quarks evolve through the emission and absorption of gluons. 
The smallness of as is never assumed, so the evolution is not perturbative. With all partons taken into 
consideration globally, we follow the evolution of the configuration in the color space. When the configuration 
exhibits color neutral subclusters, we allow branching to take place with a possible contraction of the cluster 
size in accordance to a reasonable rule consistent with confinement dynamics. Successive branching leads 
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to smaller and smaller clusters, until they are finally identified as particles and resonances. The decay of 
resonances are also taken into account before the final state of an event is determined. Evidently, the model 
contains many features of soft interaction that are familiar and desirable at a qualitative level. Here we put 
them on a quantitative basis. 

We shall call this model ECOMB, which stands for eikonal color mutation branching. Only the eikonal part 
overlaps with an earlier model, called ECCO, an eikonal cascade code |Q. Whereas branching is put in by 
hand in ECCO, it is a consequence of the color dynamics in ECOMB. Since the partons play a fundamental 
role in this model, it is significantly closer to QCD than any eikonal model on soft hadronic collisions has 
ever been. 

There are a few parameters in the model. They are adjusted to fit a large body of experimental data on 
low-pT processes for 100 GeV. They include aei, <Jinei: (n), Cq, dn/dy, and Fg for all s in the range 
10 < ^/s < 100 GeV, and all rapidity intervals. Fq are the normalized factorial moments, whose power-law 
dependence on the rapidity bin-size 6 has been referred to as the intermittency behavior |^. Except for 
Fq, all the other pieces of datga are global in nature; i.e., examined in or averaged over all rapidity space. 
They can be fitted by many models. Fq in small rapidity intervals exhibit local fluctuations, which are what 
invalidate most of those models. 

There is one important aspect about our model that should be commented on in these introductory 
remarks. ECOMB is built upon the parton model, and therefore represents an approach to multiparticle 
production at low pT that is opposite to that of the more familiar models, such as the Lund string model j^], 
and the dual parton model . The basic idea of those models is that the momentum of an incident hadron 
is carried mainly by a few of the valence quarks, which upon collision drag a color flux tube or a string 
that subsequently fragment due to confinement forces. The parton model, on the other hand, as originally 
proposed by Feynman for soft processes |pd| , regards the hadron momentum as being carried by all partons, 
whose momentum distribution is frozen by a hard collision process (except for evolution). For production 
at low pt the partons cannot be regarded as being momentarily free. It does not mean that the parton model 
is invalid; it only means that the model cannot be naively applied, as in hard processes. In describing the 
inclusive cross section of pp collision in the large xf region, the Lund model had one extreme view in 
that the produced hadrons are the result of string fragmentation, while the recombination model had 
the other extreme view in that the partons that are originally in the incident hadrons recombine to form 
the detected hadrons. Whereas the Fritiof model 0 is a refinement of the original Lund model, ECOMB 
represents an enormous step of upgrading of the original recombination model. Most importantly, it treats 
hadronization in the central region with due regard to hadron size as well as color dynamics. 

The paper is organized as follows. In Sec. 2 we review the Geometrical Branching Model (GBM) with 
emphasis on the eikonal formalism for multiparticle production. Then the parton model is incorporated into 
the GBM in Sec. 3. The parton number distribution at each impact parameter is used only as the initial 
condition for the dynamical evolution process of color mutation, which is discussed in Sec. 4. The comparison 
of the results from Monte Garlo calculation with the data on intermittency is carried out in Sec. 5. Some 
concluding remarks are made in the final section. 


II. A REVIEW OF THE GEOMETRICAL BRANCHING MODEL 

The aim of GBM is to describe multiparticle production in hadronic collisions through soft interaction 
0 - To focus on processes in which hard subprocesses are unimportant, we confine our attention to the 
energy range 10 < y/s < 100 GeV, which covers the GERN ISR energies. The model is constructed to 
possess the properties of geometrical scaling [§ and approximate Koba-Neilsen-Olesen (KNO) scaling 
that are observed in that energy range. It consists of two parts: the geometrical part that is based on the 
eikonal formalism of hadronic collisions, and the branching part that describes how hadrons are produced. 
The details of the later will be modified and improved in the following sections. Here we review the general 
framework of how the two parts are put together. 

In terms of the eikonal function H(6), which is assumed known from an independent source, the elastic, 
inelastic, and total cross sections are 


aei= I 


( 1 ) 


2 


= / cfb{l-e 


(2) 


^tot = J d^b2il-e-^^^^) . ( 3 ) 

In the energy range stated above there is geometrical scaling; i.e., Oeijotot is roughly constant That 
property can be guaranteed, if fl depends only on the scaled impact parameter i?, where 


R = b/bo{s) , (4) 

so that (0)-(i may be written as 

pOC 

(Te/= 7r&o(s) / dR^ {I - ^ (5) 

Jo 

nOO 

CTinel = ^bl{s) / dR^ (1 - , ( 6 ) 

Jo 

poo 

<Jtot = TTbl{s) 2(1 - . ( 7 ) 

Jo 

In order that the inelasticity function 

g{R) = 1 - (8) 

satisfies the normalization condition 

poo 

/ dR^g{R) = l , (9) 

^0 

we are free to set the scale bo{s) by requiring 

0■^nel = T^blis) . ( 10 ) 


The function g{R) describes the probability of having an inelastic collision at R. 

The multiplicity distribution Pn is a result of sampling the multiparticle production processes by many 
collisional events, each of which may have a different impact parameter. Hence, should have the form 

Pn = J dR^g{R)Qn{R) , (11) 

where Q„(i?) is the probability of producing n particles at impact parameter R. Both Pn and Q„ are 
normalized by 

^P„ = ^Q„ = 1 . (12) 

n n 

The eikonal formalism is, of course, quite general. It satisfies unitarity, which relates elastic and inelastic 
scattering amplitudes. Furthermore, it emphasizes the spatial properties of the colliding hadrons, which are 
known to be extended objects. Any treatment of the collision process that ignores the geometrical aspects 
of the hadrons leaves out some part of the physics of the problem, which cannot be unimportant for soft 
processes. 

The inelasticity function g{R), defined in (^), can be expanded in a power series 

OO 

g{R) = '£7r^{R) ( 13 ) 
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where 


= (14) 

/i! 

The ^th-order term may be regarded as the /rth-order rescattering contribution, and can be related to the 
^-cut-Pomeron 11,0- In this paper we use /i-cut-Pomeron only as a terminology in reference to the /rth 
term in the expansion dil)- For /r > 2, TT^ is small compared to tti except at small R for any reasonable 
fl{R). For pp collisions the well-determined form for the eikonal function is 

n(i?) =-£n (l - , (15) 

which has been used to give a good description of da/dt 0- To get an estimate of the average p, let us 
define 


p{R) — ^ (16) 

^—1 ^—1 

and 

(Ai) = J dR^g{R)fl{R) = j dR^2Vl{R). (17) 

Using (|^ in ( 0 ). one gets {p) = 1.6 Q. Thus unless i? ~ 0, high-order terms involving p> 2 are negligible. 
But for the high-n behavior of the multiplicity distribution small R collisions are very important. 

As far as the geometrical aspect of the collision problem is concerned, the above discussion summarizes 
all that is important. To go further in multiparticle production, it is necessary to model the dynamics of 
soft interaction. Before we enter into the details of that, we consider here the connection between particle 
production and geometry. At each i?, 7r^(i?) gives the probability of having /i-cut-Pomerons, for which we 
use to denote the probability of producing n particles. Thus the overall Pn is a convolution 

/ OO 

(18) 

with = 1- Going back to the general formula 0). it is clear from (18) that Q„ is a weighted average 

ofB((: 


Qn{R) = 


;^4i?)By^7rEi?), 


(19) 


fi=i 


ti=i 


where ( 0|) h as been used. 

Using (|14|) and (W), the dependence of T^/iiR) on R is shown in Fig. 1 for various values of p. Evidently, 
at large values of R only the p = 1 term is important. At R = 0, terms up to /i = 6 can still make small 
contributions. The sum in (|l3|) , when truncated at /r = 6, is shown by the heavy dashed line in Fig. 1, and 
well approximates g{R) shown by the heavy solid line in Fig. 1. Thus we see in ( p^ that at smaller R more 
terms in the sum over p are needed, resulting in more particles produced, as is physically reasonable. 


III. INCORPORATION OF THE PARTON MODEL 

In an earlier version of the GBM Q, a self-similar branching process is assumed for B((, whose n dependence 
is determined by simulation. The defect of that approach is that the dynamics of branching is ad hoc, not 
based on QCD. The guiding principle is self-similarity, i.e., a cluster of mass m branches into two subclusters 
of masses mi and m 2 with probability D(rn, mi, m 2 ) that is independent of any specific mass scale. The aim 
is to reproduce the intermittency phenomenon observed, as the resolution scale is varied |^]. We now want 
to improve the model by incorporating nonperturbative QGD to the extent feasible. Branching becomes a 
derived property in the present approach, instead of an assumed property in EGGO [Q . 
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To do so, it is necessary to put the model in a framework in which the interactions of quarks and gluons 
can be examined in detail, with hadrons being the end product of the hadronization process. Since our aim 
is to study soft interaction in low-px production processes, rigorous calculation from first principles is not 
possible. Even the use of the parton model may be questionable when the virtuality that is relevant to 
the process is not high. However, bearing in mind that Feynman’s first paper on the parton model addresses 
specifically the soft interaction problem [pd]| , we shall similarly frame our investigation in the language of 
partons without using perturbative approximations. Indeed, it is because the interactions are soft that we 
have to follow the time evolution of the partons. 

An important departure from the usual application of the parton model is that we combine the property of 
hadrons being spatially extended objects with the property of hadrons being made up of smaller constituents. 
Clearly, for most collisions whose impact parameters are nonzero, only portions of the partons in the incident 
hadrons overlap and interact, so the number of hadrons produced depends on that overlap. The event-to- 
event fluctuation of the particle multiplicity therefore depends strongly on the impact-parameter fluctuation, 
and the event averaged parton distributions, as determined from the structure functions, are of no use in 
that respect. 

To implement our amalgamation of the two properties stated above, we adopt the eikonal description of 
Sec. 2 and modify (19) to refer to parton number instead of hadron number. Thus let us define Qi,{R) to be 
the probability of involving v partons in soft interaction in the overlap region when the two incident hadrons 
are separated by a scaled impact parameter R; let be the probability of having jz partons in the /r-cut 
Pomeron. Note that we have used Qi, and in place of Q„ and Bjj, when we refer to the i/ partons instead 
of the n particles. The parton and particle distributions are, of course, very different. As in (H), Qi, and 
B^ are related by 


OO 

QAR)=g-\R)Y,^^^iR)Bf^ ■ ( 20 ) 

fl^l 

We emphasize that in (^) the subscripts v refer to the number of partons that initiate the evolution process 
leading toward the hadrons produced in the central rapidity region. The nonoverlapping portions of the 
incident hadrons produce the leading particles in the fragmentation regions. If we use • ■} to denote 

the evolution process that takes the ly partons to the n hadrons, we may express the relationship between 
Qu[R) and Q„(i?) symbolically as 


Cln{R) = £.^n{Q.{R)) ■ (21) 

The description of £u^n is the main task of this paper. 

A way of thinking about 5^ is to consider the p = \ case, for which we have a one-Pomeron exchange 
diagram for the elastic scattering amplitude with the Pomeron being cut to expose the internal lines on 
mass shell. If the Pomeron is represented by a ladder of gluons with quark loops, then cutting it reveals 
the gluons, quarks and antiquarks, which constitute the partons of 5),. Being the probability (i.e., absolute 
square of the amplitude for p + p —> v partons) of having v partons in the one cut-Pomeron, Bl is, of course, 
more general than any specific model approximating that Pomeron by ladders or otherwise. Similarly, 
is the corresponding probability, when p Pomerons in the elastic amplitude are cut. 

Since there exists no rigorous derivation of B^, we shall assume that it is Poisson distributed around 
some mean number D of partons. It is reasonable, since it is known that a cut ladder corresponds to a 
multiperipheral diagram for a i/-particle production amplitude, whose rapidity distribution is uniform, and 
multiplicity distribution Poissonian. The mean number i/(/i, s) can depend on both p and s, since we expect 
D to increase with p and s. Lacking any information on that dependence, except that the s dependence 
should be logarithmic, this being log s physics, we adopt the following parameterization 

D{p,s) = v{s)p°'^‘''^ , (22) 

where 

v{s) = + vi Ins , (23) 

a(s) = oo -|- oi Ins -I- 02 In^s . (24) 
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We shall use these in 


= , (25) 

which in turn is used in ( po|) to determine Qv{R)- 

In the next section we describe the algorithm for parton evolution in color space event by event. For each 
event we simulate a scaled impact parameter R. For that R, we use Qv{R) to simulate an initial parton 
configuration involving v partons. After the completion of the evolution process, symbolized by £v^m n- 
particles are produced for that event. Repeated simulation of many such events results in a multiplicity 
distribution at each s, a process expressed analytically by From Pn{s) we can calculate the average 
multiplicity (n) (s), and the standard moments 


C,{s) = {n^)/{nY 


(26) 


We vary the parameters in ( p^ and (|^ ) to fit the data on (n) (s) and Cq{s) for s in the range 10 < yi < 70 
GeV. 

Although the hadronization procedure has not yet been described, we give the result here first to close 
our discussion on the initial parton distribution. We do this so as not to lose sight of the global features of 
multiparticle production, when our attention is turned to the local properties in the following sections. The 
data that we want to fit are {n) and Cq from Refs. |19|. They are shown in Figs. 2 and 3. The solid lines 
are our calculated results using 


vq = —5.10 , vi = 4.03 , 

ao = -l.l, 01 = 0.41, 02 = -0.025, (27) 

in (p^ ) and (pj) with s in units of GeV^. The relative s-independence of Cq{s), which is a manifestation 
of the KNO scaling, is nontrivial when (n) (s) more than doubles over the energy range considered for soft 
production only. The fits in Figs. 2 and 3 are evidently very satisfactory. 

Before we go into the details of the color dynamics that is responsible for the calculated results, let us 
show the parton number distributions Qv{R). They are plotted in Fig. 4 for four representative values of 
R and for yfs = 52 GeV. Note that the peak does not move below iz = 25, as i? is increased; it is due to 
the fact that the minimum is 1, no matter how large R is. Note also that the typical parton numbers are 
not large by the standards of hard processes at much higher energies. That is connected with the nature of 
those partons that initiate the evolution process, a subject we turn to next. 


IV. COLOR MUTATION 

We now consider the color dynamics of soft interaction. The initial state is that there are v partons 
distributed in some fashion in a linear array in rapidity space. This is a consequence of the parton model 
in that the partons are in the incident hadrons to begin with, and the collision rearranges those partons 
and sets off the evolution process that takes those partons to the final state, where the produced hadrons 
are decoupled. The evolution process involves quarks and antiquarks emitting and absorbing gluons. Since 
the process is not perturbative, there is no analytical method to track the time development of the process. 
Thus we shall use Monte Garlo simulation to generate the configuration at each time step. 

There are two spaces in which we must track the motion of the color charges. One is the two-dimensional 
color space; the other is the one-dimensional rapidity space. The latter can be extended to include the 
azimuthal angle (j) and the transverse momentum \pt\, but in our first attempt here we integrate over those 
variables and examine the simpler problem of a ID system. As for the color space, it is 2D, and can be 
spanned by the (hjY) axes, as in SU{3) flavor. The configuration in the two spaces are coordinated in 
the following sense. Starting from the extreme left end of the rapidity (rj) space, for which the system that 
contains no partons is by definition color neutral and therefore is represented by a point at the origin of the 
color space, we move in the positive direction in rj space until we cross a color charge. At the point to the 
right of that color charge, the corresponding point in the color space jumps from the origin to a position that 
represents the color of the charge that has been crossed. Let us denote that jump by a vector. Then each 
time we cross a color charge in the rj space, there is a corresponding vector added to the previous point in 
the color space. The succession of additions of those vectors, each one starting from the tip of the previous 
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one, forms a path. Since the whole parton system is color neutral, by the time we have moved to the extreme 
right in the ij space, the path in the color space returns to the origin, thus forming a closed loop. Such a 
loop may be self-intersecting at various points. We call such a closed path a configuration of the system in 
color space. 

The essence of our evolution process is to track the configuration of the system at each time step. The 
configuration changes because the colors of the partons can mutate, as they interact through the emission 
and absorption of gluons. The way that each color charge mutates is determined by the use of an energy 
principle, which is proposed on the grounds that the system attempts to lower its energy by changing its 
color configuration, but fluctuating forces at the microscopic level cause the changes to be erratic. At 
the same time, the confinement force tends to reduce the distances among the color charges, which are 
separated initially because of the initial momenta that the partons possess, while being constituents of the 
incident hadrons. Since we work in the ID t] space and cannot use perturbative QCD, we do not treat the 
scattering problem among the partons in either the momentum or the coordinate space. The evolution of 
the configuration in the color space is what is important. 

When the closed path in the color space evolves to a configuration in which the path crosses the origin, a 
point in the evolution is then reached that requires special attention. At that point, the configuration may 
be regarded as having two closed paths, corresponding to two color neutral subsystems. We assume that a 
subsystem with no net color charge (i.e. I 3 and T = 0 in color space) is in the lowest SU{3) representation, 
that is, a singlet. This assumption is based mainly on the expectation that a nonsinglet system requires 
more energy to sustain it than a singlet system. Relaxation to the lowest energy state, despite fluctuations, 
is the principle that guides the evolution. Since two spatially separate singlet subsystems do not interact, 
we regard a branching process to have taken place: the original large cluster of partons has partitioned itself 
into two smaller color neutral subclusters. Thereafter, we repeat our procedure of color mutation, but now 
applied to the two subclusters separately. This process is repeated again and again until all subsystems 
consist of only quark-antiquark pairs that cannot branch any further. However, a qq subsystem can be a 
resonance, in which case it can decay into pions and kaons. 

So far, we have described only qualitatively the evolution process of color mutation and branching. The 
quantitative implementation of the procedure will be given below. The conceptual basis of our approach 
to soft production incorporates many features of strong interaction that are generally regarded as being 
reasonable and physical, but lacking calculational schemes with results for comparison with data. We now 
amend that defect. The procedure that we adopt to carry out the calculations is not unique, and can be 
modified and improved in time. But the framework of ECOMB is fixed by the combination of eikonalism, 
color mutation, and branching. Our present procedure is as follows. 


A. Initial Distributions of the Partons 

As described in Sec.3, we start out with v partons for a given event. It is a number generated by use of 
the eikonal formalism with variations in the impact parameter taken into account. How the v partons are 
distributed in rapidity is drawn from the parton model. An important basis of our model is that partons are 
in the incident hadrons before the collision, and are not to be regarded as qq pairs that are excited from the 
vacuum due to the stretching of strings in the string model. Thus our partons have a momentum-fraction 
distribution as determined in leptoproduction, but with as small as can be achieved. 

In the parton model there has always been a question of how gluons hadronize in low-px processes, when 
the produced particles are predominantly pions. In the recombination model for particle production in the 
fragmentization region II, the problem is handled by saturating the sea, i.e., by converting all gluons to 
quarks and antiquarks at the outset. The normalization of the hadronic inclusive cross section turns out to 
be correct, when the recombination mechanism is applied to all quarks and antiquarks that are nearby in 
phase space. We now use the same saturation of the sea in our problem here, as we attempt to treat the 
central region. Thus when we refer to partons in the following, we shall mean only quarks and antiquarks. 
The V partons are therefore color triplets and antitriplets of equal numbers so that collectively they form a 
color neutral system. Their distributions in rapidity and color spaces are to be considered separately. 
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(1) Distribution in rapidity space 


By rapidity we mean space-time rapidity ry, where 



t + z 


t — z 


(28) 


Momentum rapidity y will not be used throughout the calculation until the end, when the momentum of a 
produced particle is to be determined. We use the rj variable because the spatial separations between partons 
change during the evolution process due to the color forces. 

The usual parton distribution in momentum fraction (at the lowest accessible implies roughly a flat rj 
distribution with rapid damping toward zero at large rj. We approximate that distribution for the forward¬ 
going incident hadron by 

p(f')(r])=po, 0<V<Vc, 

Po - . Vc < fl < Vuiax , (29) 

Vmax Pc 


where po = izfipmax + Pc), which corresponds to a total of v 12 partons when integrated over rj > 0. For the 
partons belonging to the backward-going incident hadron, the ry-distribution is 


pW(7y)=po, -Pc<P<0, 

Praax P , , 

Po - , -Praax <P < -Pc ■ 

VmsLX Vc 


(30) 


The partons in each of and p^^^ should add up to be color neutral. The overall parton distribution in 
the initial state is then 


pip) = p*'^\p) +p^^Hp)- 


(31) 


For each event, the partons are, to first approximation, distributed randomly in ry according to p(?y). 

The above distribution describes the rapidities of the partons before the collision of the incident hadrons. 
Upon collision those partons interact, even if only softly. The conventional wisdom about soft interaction 
is that the range of interaction in rapidity is short. That follows from the known fact that the correlation 
between produced particles is short-ranged. Thus apart from a limited range (—ryo < P < Po) in the 
central region where the partons from the two incident hadrons interact, most partons outside that range are 
undisturbed by the collision. We describe the effect of the interaction by modifying the initial distributions 


pi'^i(p) and pi^^(p) in that interaction region to p^^^(p) and p^^i(p), respectively, where 

p(i)(p) = 6»(p -k po) 6»(po -p)+ P^^Hp) S{p - po), (32) 

2po 

p(2)(p) = Po ^(P + Po) 6>(po - p) + P^^^(p) 0{-p - Po). (33) 

2po 

Note that the total parton distribution is unchanged, i.e., 

P^^Hp) + P^^Hp) = p(p) • (34) 

The shapes of p*'^^(p) and p^^^(p) are now trapezoidal: some of the partons initially between rj = 0 and po 


in p^^\p) are dragged to the negative p region, and some between ?yo and 0 in p^^\p) are dragged to the 
positive p region. Fig. 5 illustrates how the distributions look like at a low energy where the trapezoids 
are almost triangular. This rearrangement may not seem significant for the inclusive cross section of the 
produced hadrons, since p(p) is unchanged, but it turns out to be important for multiplicity fluctuations in 
smaller bins. The point is that p^-^^ and p^^^ are two color-neutral clusters, coming from the two incident 
hadrons. The rearrangement caused by the collision creates an overlapping region so that there is no smaller 
segment within the entire range, — Pmax < P < Pmax, where one can find color neutrality, except by accident. 
The color mutation process is sensitive to the fact that the partons distributed according to p^^^(p) + p^‘^\p) 



has no color-neutral subclusters, so that the evolution must begin from the whole cluster, not from the two 
neutral subclusters, and p ^^'>, separately. 

One can find a rough parallel between our p'^^'‘^\ri) and the strings in the Fritiof model ||^ or in the dual 
parton model [Q. However, a string in DPM is stretched between a quark in one incident proton with a 
diquark in another proton. Our is a color-neutral cluster of partons distributed from —pQ to rymax- 

To summarize, whereas p^^^ {p) and p^^^ {p) are the initial parton distributions of the incident hadrons, 
p^^'^ (p) and (p) are the parton distributions after the hadrons pass through the interaction region. The 
evolution of the color dynamics begins thereafter starting with p^^\p) and p^^\p). 


(2) Distribution in color space 

In the color space spanned by and Y, a quark is represented by a vector, which has the coordinates of 
one of the triplets: (1/2,1/3), (—1/2,1/3), and (0, —2/3). An antiquark is represented by a vector directed 
opposite to one of the above. As a distribution of partons in the p space is generated, say, (p) between 
—Pmax Y p < pq^ we assign to each parton a color vector consistent with the requirement that the quarks 
and antiquarks come in pairs, but their orderings in the p and the color spaces are totally random. The 
partons in p*^^^(p) are distributed similarly, but completely independently, between —po Y p < pmax- The 
total distribution p(p) in —Pmax Y p < pmax is the sum of these two sets, whose partons are merged in the 
p and the color spaces. Thus as we move from the extreme left in the p space in the +p direction, we begin 
at the origin of the color space, and, each time when we move over a parton in the p space, we add a vector 
(corresponding to its color) with its tail to the head of the previous vector (or to the origin if it is the first 
parton). A succession of such vectors forms a trajectory in the color space. The trajectory eventually ends 
at the origin, thus forming a closed path, when p reaches the extreme right. 

Note that the path in the color space is dual to the “path” in the ID p space in the sense that the segment 
between two partons in the p space is mapped to a point in the color space, while the point where a parton 
is located is mapped to a segment of the path in the color space. Thus the closed path is made up of v 
short segments. An example of a path in color space is shown in Fig. 6. The sequence of color charges 
represented by the vectors is indicated by rbrg ■■ - br in the figure. We have used the notation r = (1/2,1/3), 
g = (—1/2,1/3), and b = (0, —2/3). The origin is labeled by 0. 


B. Evolution in the color and rapidity spaces 

We now consider the evolution of a configuration due to QCD dynamics. A nonperturbative treatment 
of ly simultaneously interacting color charges is, of course, too difficult to contemplate here. We reduce the 
problem by considering pairwise near-neighbor interaction via the exchange of a gluon in any of the s-, t-, 
or u-channel, whichever is applicable. Starting from the extreme left in the p space, we regard the ordered 
chain of v partons as having v — 1 links (with varying link lengths). Its dual path in the color space has 
ly vectors, meeting at — 1 vertices, not counting the endpoints, which are located at the origin. Pairwise 
near-neighbor interaction means that we consider the v — 1 links in p space one at a time, according to a 
rule to be specified below. After the interactions at nWv—l links are considered, the evolution of the whole 
configuration is regarded as having taken one time step, and the process is then repeated. 


(1) Color interaction 

Consider any link in the chain, i.e., a pair of partons. Suppose that the pair is (r,g). The only possible 
outcome of an interaction by the emission and absorption of a gluon, besides remaining as {r,g), is {g,r), 
since the sum of the two color vectors must remain invariant. If the pair consists of quarks, such as (r,g), 
the only outcome of an interaction is again either (r,g) or {g,r). Only in the case of a color-neutral pair, 
such as (r, f), can the outcome be any such pair, i.e., (r, f), (p,p), ib,b) or their exchanged pairings. To 
determine the specific outcome of an interaction, we use a statistical method. 

For every global configuration a of the u partons, there is an associated energy Ea- How Ea is determined 
will be described in the following subsection. Whenever a local pair of partons interact, the outcome may 
or may not affect the global configuration. Let c be the total number of possible configurations. From the 
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foregoing discussion c can only be 2 or 6, the latter being for a color-neutral pair. Our statistical approach 
to the determination of the global configuration a, consistent with favoring the lowest energy state, requires 
that the probability for configuration a to occur is 

Pa = /Z, Z = ^ , (35) 

a—1 


where /3 is a free parameter. 

Using the Metropolis algorithm, we use Pa to determine the outcome of a local interaction at every link. 
The result affects only one vertex of the closed path in the color space, one vertex at a time. As we go 
down the chain, every vertex is then tested for possible changes. At the end of one complete sweep for one 
time step, the whole closed path has undergone a color mutation, resulting in a new path. Note that our 
procedure incorporates several features of the dynamics governing the system. Locally, there is the QCD 
dynamics of gluon exchange between quarks and antiquarks. The outcome of the local interaction depends 
on the global configuration in the color space. The statistical treatment takes into account the fluctuating 
nature of the many-body problem. 


(2) Energy of a color configuration 

The color mutation process is an attempt by the dynamical system to relax to its lowest energy state. It 
starts initially at a high energy state because the collision puts the partons in a spatially spread-out domain 
by virtue of the initial momenta of the partons that are essentially uninterrupted by the soft collision, i.e., 
no large momentum transfers. Color charges that are dispersed in space have high potential energy. Since 
soft interactions do not result in significant changes in the parton momenta, the energy of the system is 
controlled mainly by the potential energies residing at all the links in the z/-parton chain. 

Specifically, consider the ith link in a chain, where i = l,2,'-',i/—1. The net color charge that the link 
sees on the left side is 


Ci — Cj , (36) 

where Cj is the vector in color space representing the color of the jth parton. Due to the color neutrality of 
the whole system, the net color that the link sees on the right side is —Ci. Since the potential energy at the 
link i is proportional to the net color charges on the two sides of the link, we write the energy Ei of the ith 
link in the simple form 


E,= 



(37) 


where the proportionality factor is absent for the following reason. The confinement potential should be 
proportional to the distance between color charges. Here Ci is the total color of all the partons to the left of 
the Ah link, so the relevant effective distance should be measured from the center of those partons on the 
left side to the center of the partons on the right side. That is a global distance that is roughly half the total 
rj range, insensitive to the position of the Ah link. Since Ea is never used without /3 in (p5|), the omission of 
the proportionality factor in (H) is equivalent to letting it be absorbed in the definition of /3. The energy 
associated with the whole configuration a is 


v-\ 

Ea=Y,E, , (38) 

i=l 

where the sum obviously depends on the particular path in color space that the configuration a takes. 

It is this energy Ea in ( |38|) that the system attempts to lower by color mutation, as it changes from one 
configuration a to another, a'. Because of the additive nature of Ea in (^^, the probability of mutation is 
determined by the local energy difference, as can be seen in the following example. Consider a link for which 
the number of possible configurations is two, i.e., c = 2 in (|^). Then Pa can be written as 
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1 


(39) 


Pa = 


where AE is the change of energy in going to the new configuration o' 

AE = Ec' - Pa = E[ - P,. (40) 

P' is the energy at the ith link in the new configuration. Note that AE is now locally 
link i. However, the energy Ei and P' depend on the global configurations in the color 
and d^). 

Since the net color charge of all the partons on one side of a link can be large, P^ 
according to (|3^). Thus unless /3 is a small number, (|3^) can give a probability distribution Pa that can 
have an abrupt change at AE = 0. To round off the edges, /3 would have to be quite small. 


determined at the 
space, due to (^) 


can also be large. 


(3) Spatial fluctuations 

The length of the zth link in the rj space is 

di = r]i+i - rj, , ( 41 ) 

where rji is the rapidity of the ith parton, counting from the left at 77 = —rjmax- For a spatially extended 
color-neutral system, we have positive link lengths, dt > 0. We put dt = 0, when two partons are in the 
same rj bin. The rapidity range, —jymax < ^ < ? 7 max, is divided into 256 bins in our numerical calculation, so 
the smallest bin size is 277r[iax/256, which, of course, depends on the collision energy. As far as the problem 
regarding the spatial distribution is concerned, we shall assume energy independence and work only in the 
space that has 256 basic unit^whatever r^max is- We denote the width of the basic unit by Sr;. 

On the basis of (^) and (|^ ) the system undergoes color mutation, as it evolves toward the final state. 
In addition to the fluctuation of the color charge Ci, the length of the ith link, di, which is an integer in 
units of 6r], can also fluctuate. Whether the link length contracts or expands depends on the attractiveness 
or repulsiveness of the net color forces that act on the two ends of the link. To determine the nature of 
that force is beyond the scope of this treatment. We shall model the change in link length by a stochastic 
approach, consistent with how we have handled the color mutation part of the dynamics of the complex 
system. We allow di to change by an integer m^, i.e., di ^ di + rrii, where the probability for rrii is specified 
by a distribution Pi{m). Apart from the consideration that a contraction cannot render di negative, we 
would set Pi(rn) to be a uniform distribution from —rri- to to+, where m± > 0. However , if to_ > di, 
we truncate Pi{m) at m = —di and place the probabilities for contraction from —di to —rri- at m = —di. 
Thus our two-parameter formula for Pi{m) can be rewritten in the form that has a constant term plus a 
Kronecker delta term at —di. That is, with 71 and 72 as the two revised parameters, we have 

Tiim) =-fi0{m + di)9{'j2 - m) + [1 - ^ 1(^2 + di)] 6m -di , (42) 

which satisfies = 1- This distribution allows m to be positive (expansion) uniformly up to 

72 , and negative (contraction) down to —di. The net probability for expansion is 7172 , and for contraction 
is 1 - 7172 . 

Confinement implies that there should be a net contraction between the subcluster of charge Ci and 
the opposite subcluster of charge —Ci. However, these two subclusters are at an effective distance much 
larger than di, which is the length of the link specifying the separation between the boundaries of the two 
subclusters. At the ith link, there are microscopic color forces that can be attractive or repulsive. Vi{m) is 
a two-parameter description of that stochastic force, which should lead to a net macroscopic contraction for 
the whole cluster. That property is related to the value 7172 , which should be less than 0.5, if contraction 
dominates over expansion. Of course, this change of di should be tested at all links in the neutral cluster 
before a net effect is known. 


C. Branching 

As we have discussed in the introduction of this section, the evolution process described here is a branching 
process. A fission of the color-singlet cluster occurs, when it contains two color-singlet subclusters, since 
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no confining force exists between them. The partitioning takes place when the closed path of the mother 
cluster evolves to a configuration where the path passes through the origin in the color space so that there 
are two color neutral subsystems. Thereafter, the color mutation process is applied to the two subsystems 
separately and independently. This is repeated again and again until all subclusters consist of only qq pairs. 

During color mutation, it is possible that a baryon-like singlet can be formed with a (r, g, b) closed loop. 
Since baryon-antibaryon pair production in pp collision is rare at energies where soft interaction is dominant, 
we shall ignore 3q and 3q singlets, (since mass inhibition is not explicitly taken into account), and require 
color mutation to continue until only mesons are produced. If in the middle of a color chain there is a color 
neutral qq pair, such as rf, it is not allowed to break off from the chain; it must continue to mutate until a 
genuine neutral subcluster, counting from the extreme left or right in the iq space, is formed. 

In Fig. 7 we show an example of how 24 partons evolve under color mutation, and how branching takes 
place. Thirty-two time steps are taken to get to the final state of qq subclusters only, when hadronization 
occurs. 

To gain some insight on how long a typical evolution process takes for a cluster to partition into two, we 
show in Fig. 8 (a) the distribution in time steps for 10^ initial clusters each consisting of ten pairs of qq with 
random color ordering. The average is about 3.5 steps for just one branching. For the evolution processes 
to complete, the daughter clusters must continue to branch successively, and the overall distribution at the 
end of the evolution is shown in Fig. 8 (b). Note that on the average quite a large number of time steps is 
required for the evolution of an average-sized cluster, much larger than what is usually involved in a hard 
porcess where high virtuality is degraded by branching. 

At every stage of the evolution process, a cluster shrinks due to overall spatial contraction. We always 
keep the center of the cluster invariant to conserve momentum. When a branching occurs, the two daughter 
clusters will have their own respective centers, which will remain invariant during contraction, until they 
themselves branch. At the end when a hadron is formed from a qq pair, the hadron momentum rapidity y 
will be identified with a value of g taken randomly between the 77 values of the quark and antiquark. The 
reason for doing this is that without knowing the mass and the transverse momentum px of the hadron, a 
precise determination of its y is not possible, nor meaningful. A random value within a small range is good 
enough. A rough identification of y with g is justified for free particles at high energy. We do not expect 
that our measure of fluctuations in the final result will be sensitive to this y-g identification. 


D. Hadronization 


The branching process terminates, when all the neutral subclusters are reduced to the composition of qq 
pairs only. Those subclusters are identified as hadrons, which may be pions, kaons, or resonances. Those 
resonances must be allowed to decay before the total number and distribution of particles are counted for the 
final state of the event. The probabilities of producing various resonances and stable (in strong interaction) 
particles have been studied experimentally in pp collisions in [ p0[ . We use that reference as a generic guide 
for the proportions of all particles produced in any general hadronic collision. 

In our simulation we use that guide to determine whether a neutral subcluster is a resonance or not. In 
either case we give it a transverse momentum according to a specified distribution. Let us use kx to denote 
the transverse momentum of a qq pair. The distribution we use is 


dn 1 _i. 2 //j. 2 \ 

dk'^ (fcy) 


(43) 


where is a parameter to be varied. If the qq pair forms a pion or a kaon, then that kx becomes the px 
of the particle without any change. However, if the qq pair forms a resonance, then we assume an isotropic 
decay distribution in the rest frame of the resonance. The azimuthal angle (jj of k is assigned randomly. 
After boosting back to the cm system, the 3-momenta p of the decay particles are then determined. The 
collection of all the final-state particles results in the exclusive distribution in the {y,px) space for the event 
generated. 

The charges of the particles are chosen randomly subject to the constraint that the total charge can only 
be either 0 , 1 , or 2 , on the grounds that the leading particles in a pp collision can be doubly-, singly-, or 
un-charged. The particles that we treat are all the produced particles, excluding the leading particles. Once 
the charge of a specific particle or resonance is generated, the decay of the resonance is then assigned the 
usual branching fractions into different channels according to the particle data book. 
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V. RESULTS 


In the previous section we have introduced seven parameters: rjc and ? 7 max in (|^) and (|30|), ryo in (|^) and 
(^, (3 in d^), 7 i and 72 in (^), and in (|4^). They are to be varied to fit the data on inclusive distri¬ 
butions and on fluctuations of the exclusive distributions. The properties of fluctuations will be quantified 
in terms of the normalized factorial moments, Fq{5) for farious bin sizes S. 

The parameters rjci Vmax, and (k^) are essentially kinematical; they set the boundaries of the phase 
space in which the partons are placed. They do not affect the dynamics of color mutation and the spatial 
fluctuation of the clusters. The data used to determine them are the rapidity distribution dn/dy 0 and the 
transverse mometum distribution g{pT) of the produced particles in the final state. The energy range of the 
data in 0 is 22 < y/s < 63 GeV, for which hard scattering is negligible. The values of ijc and rymax depend 

on s. To fit the data, we have used the values given in Table I. The value of is chosen to be 400 

MeV. In Fig. 9 we show the rapidity distributions where the histograms are the results of our calculation, 
while the data are from |21|. There, 77 refers to pseudaorapidity, not space-time rapidity. The agreement is 
clearly satisfactory. Although in carrying out our calculation the other parameters need to be specified also, 
we have not given them in Table I because the results on dn/drj are insensitive to them. They will be given 
below. 



22.0 

23.6 

30.8 

45.2 

53.2 

63.2 

Vc 

1.76 

1.55 

1.42 

0.889 

0.762 

0.635 

'^max 

5.0 

6.6 

6.6 

6.5 

6.5 

6.5 


TABLE I. The parameters for the initial condition of partons in the evolution 


The multiplicity distribution at y/s = 22 GeV is shown in Fig. 10 for \y\ < 2.5. The calculated result 
evidently agrees very well with the data |^, especially at the high n end. At the low end, there is a slight 
discrepency, which could very well be due to the fact that we have not included diffractive production. That 
portion of the distribution will have negligible effect on our calculation of the factorial moments. We have 
not included in Fig. 10 the for narrow rapidity windows, since essentially the same physical content will 
be conveyed by the Fq{5) to be presented below. We have not considered the UAl and UA5 data at higher 
energies because those data contain significant contributions from hard subprocesses. The inclusion of such 
perturbative processes is not difficult , and should be investigated in order to extend the energy range 

where EGOMB can be applied. 

The parameters that influence the evolution process of mutation and branching are 70 > Pi 7i and 72 . 
70 specifies the overlap region of the initial color-neutral clusters, P pertains to the probability of color 
mutation, and 71 and 72 characterize the fluctuations of the link lengths. The data used to determine them 
are on factorial moments and their variants ||^. It is here that we can underline the importance of those 
data on fluctuations, without which we have no guidance on how to restrict the detail dynamics of particle 
production. Putting that in another way, in the absence of a procedure to calculate from first principles, 
any model that fails to fit the fluctuation data is missing some aspect of the basic dynamics. To the extent 
of our awareness, very few models on soft interaction have been put to the test of confronting those data on 
the factorial moments for varying bin sizes. 

The normalized factorial moments, first suggested by Bialas and Peschanski can be expressed as 


. _ 1 (n(n-l)...(n-g-H))fe 

’ {n)l 


where n is the multiplicity in a bin of size <5, the average is performed over all events for the /cth bin, and 
the summation is over all bins of the same 6, M being the total number of bins {2Y/5). Data on Fq for soft 
interaction have been obtained by NA22 |^|. In recent years Bose-Einstein correlation between identical 
particles has been recognized as having an important effect on the intermittency data at very small S |^]. 
However, at this stage of our development of EGOMB, we have not included BE correlation, since it is 
unrelated to the color dynamics that we attempt to tune. We therefore limit our scope in this first attempt 
and focus here on Fq, for which we fit the data of ||^ as functions of bin width, rather than on correlation 
functions at very small momentum-difference squared. 

The parameters that we have varied to give the best fit of the intermittency data at -^/s = 22 GeV are: 
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(45) 

(46) 

(47) 


i){s) = 9.1, o(s) = 0.63, 

??max = 5, r]c = 3.5, 770 = 1.9, 

/3 = 0.0015, 71 = 0.077, 72 = 5. 

Note that 7172 = 0.385, implying that it is roughly twice as likely for a link to contract than to expand. To 
explain why the hrst four parameters listed above are different from the ones used previously for the ISR 
data (in Sec. 3 and Table I), we remark that we have not been able to derive from the values of Fq given 
by the data of NA22 the values of Cq that can agree with those of ISR . Thus while the parameters 
in Table 1 are adequate for the ISR data that range up to 63 GeV, we must use a slightly different set of 
parameters at ^/s = 22 GeV for the NA22 data that provide the only information on intermittency. 

In Fig. 11 we show the intermittency data of [ p^ , which to our knowledge have not been reproduced 
by any model, except ECCO The lines shown in Fig. 11 are the result of our present calculation in 

EGOMB. Evidently, the agreement between our results and the data are satisfactory for all values of 6 y and 
q. It should be noted that to achieve the fits attained is highly nontrivial. If any part of the dynamical 
process in generating the hadrons is altered, one would not be able to obtain the rising factorial moments, no 
matter how many parameters are used. By working with the many parts of our model, all of which affect the 
determination of Fq , we have gained conhdence in regarding the dynamics of color mutation and branching 
as having captured the essential properties of soft interaction. 


VI. CONCLUSION 

It is rather satisfying that we have been able to reproduce the intermittency data in Fig. 11. Scaling 
behavior of that type implies self-similarity in the dynamics of particle production. The branching process of 
color-neutral clusters, successively partitioning into small clusters, is self-similar, since the same algorithm is 
used at every step. Thus from the outset ECOMB has a chance of generating a scaling behavior, in contrast 
to the string model or dual parton model, which have been shown by NA22 not to fit the intermittency 
data 1^ . However, to have self-similar dynamics does not guarantee that the observed scaling behavior can 
be obtained. An appropriate balance of contraction and expansion of the clusters during the evolution is 
necessary to capture the essence of the color dynamics of a complex system. Furthermore, the initial state 
of the parton system at the start of the evolution process is also important, since it can dictate the bin 
size from which the scaling behavior commences. These physical aspects of the model must be incorporated 
properly in order to produce the result shown in Fig. 11. 

We have not considered the intermittency data of UAl because they contain minijets. If one is willing 
to adopt the point of view that such semihard subprocesses are not dominant enough to alter the nature 
of intermittency, one could attempt to apply ECOMB to the UAl data. However, it is well known that 
geometrical scaling and KNO scaling are violated at the SPS energies, so some modification is necessary. 

The most important omission in this paper is Bose-Einstein correlation, without which our model is not 
completely realistic. That defect must be corrected in an improved version. At this stage we can only offer 
this version as a first step in the right direction. Whereas the details of our parametrization here may alter 
after the BE correlation is incorporated, the branching dynamics of color mutation will hold as the basic 
framework of the model. Since the power-law behavior of Fq has been observed for nonidentical particles 
, we know that self-similarity is necessary in the absence of BE correlation. After the symmetrization of 
identical-particle states is performed, the strength of intermittency may be enhanced. The parameters may 
therefore have to be readjusted, but the main features of the dynamics need no modihcation. 

We have amalgamated many concepts that form various elements of the conventional wisdom about soft 
interaction. They include: (a) hadrons having sizes, (b) eikonalism, (c) parton model, (d) interaction of 
quarks via gluons, (e) statistical properties of a many-body system, (f) spatial contraction and expansion 
of a color system, and (g) resonance production. They are interlaced by intricate connections described 
in this paper. One may regard the approximation of nonperturbative QGD by color mutation as a gross 
simplification. We conjecture that the final realistic result will be insensitive to the details of how well many- 
body QGD is approximated, since the model is constrained by so many other aspects of soft interaction. This 
work demonstrates the importance of how each and every one of those separate components of the overall 
structure must work together in order to achieve the scaling behavior of local fluctuations that is observed. 
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# of partons : 24 


B : Blue 
b : Anti-blue 


G: Green 
g : Anti-green 


R: Red 
r : Anti-red 


1 BGGBgbBgGggbbbggGGRBBGrb 

2 BGGgBBgGggbbbbggGGRBBGbr 

3 GGBgBBgggbGbbbgGgGBRGbrB 

4 GBgBGgBgbGgbbgbGgBRGGbBr 

5 GBgBgBGbggbbGbGggBRGbGBr 

6 BGggBGbBgbbGgbggBRGbGBrG 


7 GBgBgbBGgbGb 

8 BgBgbBGgGGbb 

9 gBgbBBGgGbGb 

10 BgbBBggGbGGb 

11 BgbBBgGbgGGb 

12 BbBBggGbgGGb 


13 bB BBgGbggGbG 

14 bB BBGbgggGGb 


15 Bb BGbgBggGbG 


16 

Bb 

BbGg 

1 1 

gBGbgG 

1 1 

17 

bB 

Bb 

gG 

BGgb 

Gg 

18 

Bb 

bB 

Gg 

GBgb 

gG 

19 

bB 

bB 

Gg 

BGgb 

gG 

20 

bB 

bB 

gG 

GBbg 

gG 

21 

Bb 

bB 

Gg 

GbBg 

gG 

22 

Bb 

Bb 

gG 

GBgb 

Gg 

23 

bB 

bB 

gG 

BGgb 

gG 

24 

Bb 

Bb 

Gg 

BGgb 

gG 

25 

Bb 

Bb 

Gg 

GBbg 

Gg 

26 

Bb 

bB 

gG 

GBbg 

Gg 

27 

Bb 

Bb 

gG 

BGbg 

gG 

28 

Bb 

bB 

Gg 

GbBg 

gG 

29 

Bb 

bB 

gG 

GBgb 

gG 

30 

Bb 

Bb 

gG 

BgbG 

Gg 

31 

bB 

Bb 

gG 

BbgG 

gG 


I ^22 

32 bB Bb gG Bb gG gG 


gbggRBGGbBGr 

ggbgRGBbGGrB 

gbggRGBGbrGB 

gggbGRGbrGBB 

gggbRGGbGBBr 

ggbgGRGbBBrG 

gbgGgGRbBrGB 

bggGGRgBrGBb 


bgGGggRrGB bB 

gbGGggrGBR Bb 


gbGGgrgGBR bB 

gbGgGrgBGR bB 

gGgbGgrBRG Bb 


Gg 

gGgrBRGb 

1-1 

Bb 

gG 

gG 

grBRGb 

bB 

Gg 

Gg 

rBgGbR 

bB 

Gg 

Gg 

BrgbGR 

bB 

gG 

gG 

rgBbRG 

bB 

Gg 

Gg 

rgbBRG 

Bb 

Gg 

Gg 

rbgBGR 

Bb 

gG 

gG 

brBgRG 

bB 

gG 

gG 

brBgRG 

Bb 

Gg 

gG 

rbBgRG 

bB 

Gg 

gG 

rBbRgG 

1-1 

Bb 

Gg 

gG 

BbRr Gg 

Bb 


Gg gG Bb Rr Gg Bb 



Probability 


FIG. 7. An example of a branching process that starts with 24 partons. The numbers in the hrst column denote 
the nnmber of time steps taken for the conhguration in any particular row. 
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FIG. 8. (a) The distribution of times steps taken for 20 initial partons before a hssion process occnrs. (b) The 
distribution of time steps taken for the 20 partons to evolve to the end when branching terminates. 
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FIG. 11. Intermittency data of normalized factorial moments for g = 2 — 5. The data are from Ref. [25]. The lines 
are determined from ECOMB. 



